Narrowband Biphoton Generation due to Long-Lived Coherent Population 

Oscillations 
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We study the generation of paired photons due to the effect of four-wave mixing in an ensemble 
of pumped two-level systems that decay via an intermediate metastable state. The slow population 
relaxation of the metastable state creates long-lived coherent population oscillations, leading to 
narrowband nonlinear response of the medium which determines the spectral width of the biphotons. 
In addition, the biphotons are antibunched, with antibunching period determined by the dephasing 
time. During this period, damped oscillations of the biphoton wavefunction occur if the pump 
detuning is non-zero. 

PACS numbers: 42.50.Ar, 34.80.Pa 



Traditionally, paired photons are produced from spon- 
taneous parametric down conversion in nonlinear crys- 
tals. The bandwidth of such biphotons is very broad and 
typically in the terahertz range [l[, which makes them 
useless for some applications in quantum information sci- 
ence which require strong interaction between photons 
and atomic systems. This problem can be overcome 
by generating biphotons in cold atomic systems which 
have a narrowband nonlinear response. For example, 
biphotons can be produced when a double-A system [2- 
4] or two-level system (TLS) [g, [f| is pumped by two 
counter-propagating laser fields. Then phase-matched 
and energy-time entangled photon pairs are produced 
due to the effect of four- wave mixing (FWM) . Biphotons 
from a such source have a bandwidth in the megahertz 
range and coherence time of hundreds of nanoseconds. 

In this paper, we demonstrate that narrowband bipho- 
tons can also be produced due to the effect of long-lived 
coherent population oscillations (CPOs) in a TLS with 
an intermediate metastable state. In such a system, the 
width of the nonlinear response is determined by the life- 
time of the metastable state, which can vary significantly 
depending on the nature of the quantum system. For ex- 
ample, in semiconductor quantum wells and dots [7| the 
CPO lifetime is in the microsecond range, whereas in a 
ruby crystal Q or organic film Q it can be more than a 
millisecond, leading to a broad range of potential appli- 
cations of photon pairs based on the CPO effect. 

We consider the interaction of an ensemble of TLSs 
that decay via a single intermediate metastable state 
with two counterpropagating pump fields with amplitude 
E Q (see FigfT|. The medium is assumed to be optically 
thin in the direction of pump propagation and the ef- 
fect of pump depletion is not taken into account. Due 
to pumping of the TLS by two counter-propagating laser 
fields, photon pairs are produced [5|, |6| where photons 
from the same pair also counter-propagate so that the 
phase-matching condition of FWM is satisfied [10| [see 
Fig. da)] (actually, in such a configuration biphotons are 
emitted into the whole Air space). In order to allow for 
the spontaneous initiation process, the generated weak 
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FIG. 1: (color online). In the presence of the counter- 
propagating pump fields, phase-matched counter-propagating 
biphotons are generated inside the medium due to FWM. 



fields are described by quantum-mechanical operators 

4$ = E i%,2 (*, *) and E{; 2 } = E[%1 2 (z, t) , (1) 

where the subscript 1 denotes the field at frequency 
Ui = ujq + 8 propagating along the positive z-axis, and 
the subscript 2 denotes the field at frequency 0J2 = ujq — 5 
moving along the negative z-axis, E[ 2 — (hioi^/^EoV) 
is the vacuum field, V is the quantization volume, and a 
and a^ are the photon annihilation and creation opera- 
tors. 

To describe the evolution of the atomic ensemble, we 
begin with the Heisenberg operator equations of the mo- 
tion in the dipole approximation: 

h (d/dt + r fcQ + iu) a ba = id ba E (+) (a bb —<?„„), 

h (d/dt + T ba - iu) o ab = -idabE^ (abb -a aa ), (2) 

h (d/dt + jb) a bb = id ab E ( ^a ba - idb a E^a a b, 

(d/dt + 7 ca ) a cc = 7h c 5 fc6 , 

where i?^ = \i) (j\ is the atomic operator, E&j is the 
total field operator, and dij is the transition dipole matrix 



element, T ba is the transverse relaxation rate, 7^ is the 
longitudinal decay rate from the state \i) to the state \j), 
lb = Ibc + Iba is the total decay rate from the excited 
level, and Q — ujq — Wf, a is the pump detuning from the 
resonance. We also assume that the system is closed so 

that Oaa + &bb + <7 CC = 1 . 

We apply the slowly-varying envelope approximation 
and write the total held operator as 



£(±) 



E (.±) e TiSt , £(±) e ±«A e T^ot 



(3) 



To eliminate the fast oscillating term in Eqs. (J5J, we in- 
troduce the transformations ^ba,ab (i) = Oba,ab (t) e^™ ' 
and crjj (t) = a.jj (t) \j =a .b,c- In order to find the medium 
response to the weak generated fields, we apply the Flo- 
quet theory [ll| and write 



r(°) 
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. a (-s) e ™ 



(4) 



The zeroth-order solution of Eqs. © - (@J gives the re- 
sponse of the medium to the pump field and the popu- 
lation distribution between the quantum states, whereas 
the first-order solution determines the medium response 

to the weak generated fields since P\ t 2 = Ndba&ba 0- 
The pump-probe interaction with a TLS is character- 
ized by population beating or coherent population os- 
cillations (CPOs) at S, the frequency difference between 



pump and probe fields (12J. In an ordinary TLS, the 
CPOs decay at the same rate as the excited state. How- 
ever, the situation can be quite different if an interme- 
diate metastable state is included [see Fig. QJb)]. In 
the case where the rate of transfer of CPO from the ex- 
cited state to the metastable state is much faster than the 
rate of transfer of CPO from the excited to the ground 
state due to population relaxation or pump-induced tran- 
sitions, that is, 



Jbc > Vo,76a, 



(5) 



where Vq = db a Eo/ (2/t) is the pump Rabi frequency 
which is assumed to be real, long-lived CPOs of ground 
and metastable states are created 13J. This leads to a 



narrow dip in the probe absorption spectrum and a nar- 
row peak in the FWM spectrum [13l4l5l | . 

Under these conditions and taking into account that 
level |c) is the metastable level j ca <C jba, Ibc the steady- 
state response of the medium to the generated fields is 
given by [lj| Il6| 



J+s) 

J ba 

T (-s) 
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= (aiE[ 

= (a 2 E, 



(+) 



(-) 



+ PM ] )d ba /K (6a) 

+ /? 2 £ 1 +) ) d ab /h, (6b) 



where a±.2 are proportional to the effective linear suscep- 
tibilities and (01,2 are proportional to the effective third- 
order nonlinear susceptibilities and are responsible for 
the generation of the paired photons. They are given by 



ai, 2 = ±t (1 + -X") / [T li2 (1 ■ 

1)2 = ±tx/ [r lj2 (i + k)] , 



0] 



(7a) 
(7b) 



where 

X = -n lca / (W - iS) 
is the coherent field interaction term and 

W = {l + n) lca 



(8b) 



determines the characteristic width of the window in 
which coherent interaction between the fields occurs, 
k = 2Vq/ [jca^ba (1 + ^ 2 /rL)] * s * ne saturation param- 
eter, and Ti t 2 = Tt a + i (T^ — S) . 

The evolution of the annihilation and creation oper- 
ators a and a' is described by the coupled propagation 
equations (J): 



d d 
dt dz 
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dt' 


dz 



\a 1 (z,t)=ig 1 Na { + 5 \ 
a\(z,t) = -ig 2 Na { ab d \ 



(9a) 



(9b) 



where g\ t i = db a E\ l/H are the coupling constants. 

In order to solve these coupled equations we make a 
Fourier transformation, neglect the term iui/c as it does 
not affect the final result j4|, and substitute Eqs. (|6aj) 
and (|6b|) into Eqs. (|9"a|) and (|9bl) to obtain 

(fOa) 

(10b) 

where ( = er cs r ba N/2, a cs = dl a uj Q / (ce Q hT ba ) is the 
atomic cross section, N = N/V is the atomic density 



ddi/d (C,z) = iocidx + i0id 2 , 
da 2 /d ((z) — ~ia 2 a 2 — i/^Si, 



(0) 



p(0) 



and it is also assumed that E{ ' w E 2 ■ The biphotons 
counter-propagate so that photon '1 ' leaves the medium 
at the point z — L and photon '2' leaves at z = 0. The 
boundary conditions derive from the vacuum field fluc- 
tuations at z = for photon '1' and at z — L for photon 
'2' [see Fig. QJa)]. Thus, the solution of this system for 
variables a\ (L) and a 2 (0) of the backward-wave prob- 
lem can be written as a linear combination of the initial 
boundary values 



Si (L) = AiSi (0) + BiS 2 (L) , 

a 2 (0) = A 2 Z\ (L) + B 2 a x (0) , 



where 

A x . 2 



,±i(a 
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Bi, 2 = i/3i, 2 sin [RLj / (RD) 
D = cos RL — i — — — sin 



(Ha) 
(Hb) 

(12a) 
(12b) 
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in (rl\ , 



1/2 



R= (ai + 012) /4-0102 saidL^CL. 

The correlation between the photons emitted to the 
left and right is described by the second-order Glauber 
correlation function G 21 |17| 

G 2 2 ?(T) = (a\(L,t)al(0,t + T)a 2 (0,t + T)a 1 {L,t)y 

(13) 



As a field emitted by many statistically independent 
atoms behaves as a Gaussian random variable, we can 
use the Gaussian momentum theorem [4|, LL7J and rewrite 
Eq. (H3]) in the form 



G 



(2) _ 
21 — 



G«(0)GW(0) + |$ 21 (r)r 



where the terms 



G ! S(0) = (aI !2 (0)a 1 , 2 (0) 



(14a) 



(14b) 



describe the appearance of uncorrelated photons which 
produce a flat background, and the second term 



*ai(T) = <3a(0,t + T)oi(L,t)) 



(14c) 



describes the appearance of entangled photon pairs and 
corresponds to the biphoton wave function 0, [13, LL8| ■ 

As seeding fields Ex 2 are absent, the vacuum field fluc- 
tuations determine the initial conditions and taking into 
account the commutation relation for the input field op- 



erators 



a 1 ^ 2 (z,uj),a\ 2 (z,-uj / ) 
in Eqs. (|14bl) and (|14cl) . we obtain 
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2ttc 



L/ (2ttc) S(u} + u') 
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D 
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* 21 (T) = 2^ 
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(15a) 

(15b) 



In particular, we are interested in the biphoton coher- 
ence time which is determined by the width of the func- 
tion $21 (t). Harris and coworkers Q have pointed out 
that a long coherence time can be obtained due to the 
effect of slow light [19| experienced by one of the pho- 
tons of the entangled pair but not by the other. Here we 
demonstrate that a long coherence time can be obtained 
even in a optically thin medium 



ai 2 £, j0i,2-£ < 1, 



(16) 



where the time delay between the photons due to the 
slow light effect is negligible. 

Under the conditions of Eq. (flT)J). Eqs. ([T2"al) and (fT2b|) 
simplify to 



A 12 = 1, B 1:2 = iftipL. 



(17) 



To find the analytical form of the correlation function, 
we substitute Eqs. (JTTJ) into Eqs. (|15ap and (|15bl) and 
integrate over 5: 
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r/<5 



(18b) 




FIG. 2: (color online) Normalized second-order correlation 
function for k = 1 and 7ca/r(, a = 10~ 4 , for il = (red solid 
line) and fi/rf, a = 10 (blue dotted line). The normalization 
constant £ = |/3<5=o| L 2 has the same value for both cases. 
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and also we as- 



where the constant n = ^ , ; 77- 

sume that r& a ^> W. The normalized second-order cor- 
relation function is then given by 



$ (r) = 



G% (r) 



G\ 1] (0) G 2 l) (0) 



-W\r 



1 

1 + e 



-2r ba |r| _ 2 cos (fir) e" 



IA 



.5=0 



! L2 



(19) 



When Eq. (|16l) holds, the denominator of the second 
term of Eq. (I19[) is much less than unity and as a result 

the visibility Kis = (ffmax - ff^L) / (ffmax + ffmin) ~ L 
In Fig. @, we show the behavior of the normalized 
second-order correlation function. It can be seen that the 
characteristic width is determined by 1/W and also that 
there is an antibunching-like effect with a characteristic 
time 1/Tba (red solid line). When the pump detuning is 
non-zero (blue dotted line), damped oscillations at fre- 
quency $7 with a decay rate 1/Tb a are produced. 

This behavior is easy to understand if we look at Eq. 
(|18b[) . which shows that the biphoton wave function is a 
coherent superposition of two different FWM processes 
Q . The first term describes the generation of the bipho- 
tons with a spectral width W centered at the point S = 
(we call this FWMi) and the second term describes 
generation of the biphotons with a spectral width Tb a 
from the two sidebands which are centered at the points 
6 = ±fi (called FWM n ). As we can see from Eq. (|18bl) . 
the phase shift between FWMj and FWMjj oscillates 
at the frequency of the pump detuning which leads to 
constructive or destructive interference between them, as 
shown in Fig. ([2]). The coherence time of the bipho- 
ton generated by FWMn is determined by 1/Tf, a , so 
that these biphotons can contribute to the coherent su- 
perposition only during this period, which causes a fast 



damping of the oscillations. When the pump detuning 
is zero, there is always destructive interference between 
the biphotons from FWMi and FWMu, leading to an 
anti-bunching dip with a width l/T^a (wc should note 
that the depth of the antibunching dip never goes to zero 
as we have a macroscopic ensemble of the quantum sys- 
tems). The long coherence time is caused by biphotons 
that originate from the FWMi process as it has a very 
narrow bandwidth [see Eq. (|8b|) ]. 

In summary, we have demonstrated that the combined 
effects of FWM and long-lived CPOs in a TLS with inter- 



mediate metastable state are able to produce narrowband 
biphotons with a long coherence time whose maximum 
value is equal to the lifetime of the metastable state. The 
biphotons' waveform and bandwidth can be controlled by 
the pump intensity. During the time 1/Tb a the biphoton 
wavefunction shows antibunching behavior. If the pump 
field is detuned, damped oscillation during this period is 
observed. 

The authors thank A. Pe'er and A. Eilam for helpful 
discussions. 
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